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Goals  of  the  Research 


This  Grant  aimed  at  an  understanding  of  the  strongly 
nonlinear  electrical  properties  of  charge  density  wave  (CDW) 
conductors. 

Some  progress  had  been  made,  in  understanding  weakly 
nonlinear  properties,  with  models  which  retained  the  spatial 
randomness  of  the  pinning  potential  (due  to  crystal 
inhomogeneities) .  The  complexity  of  a  random  potential  (when 
combined  with  the  dynamic,  many-body  and  nonlinear  nature  of  the 
problem)  meant  however  that  progress  into  the  strongly  nonlinear 
regime  was  extremely  slow,  despite  the  attentions  of  numerous 
researchers.  The  PI  adopted  what  was  at  the  time  a  somewhat 
unfashionable  approach.  It  was  decided  to  see  whether  the 
randomness  in  the  pinning  potential  might  not  in  fact  be  an 
unnecessary  complication.  It  was  decided  to  attempt  to  solve  for 
the  dynamical  properties  of  incommensurate  chains.  These  have  a 
periodic  rather  that  a  random  pinning  potential,  but  they  retain 
the  collective,  or  many-body,  aspect  as  well  as  the  nonlinearity 
of  sliding  CDW's.  Once  their  dynamic  properties  are  solved  we 
must  compare  with  experiment  to  see  whether  these  models  really 
contain  the  essence  of  CDW  conduction. 

The  results  are  summarized  in  the  next  section  (II) .  A  more 
complete  account  of  the  progress  under  the  grant  is  contained  in 
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the  publications  (Section  III) ,  with  some  ongoing  work  described 


in  the  final  section  (V)  .  In  section  IV,  I  outline  a  new 
perspective  that  has  emerged  from  the  results  of  all  these 
projects  taken  as  a  whole. 

II .  Theory  vs .  Experiment 

From  the  publications  (in  many  cases  from  the  Figures)  it  is 
seen  that  incommensurate  chains  give  a  surprisingly  good  account 
of  the  following  dozen  measurements:  both  components  of  complex 
ac  conductivities  as  functions  of  field  and  frequency,  in  both 
metallic  and  semiconducting  CDW  materials;  dc  characteristics; 
scaling  of  ac  and  dc  conductivities'  elastic  properties  - 
Young's  Modulus  and  Q-f actor  as  functions  of  voltage;  bulk 
oscillations;  and  both  amplitude  and  phase  of  both  the  second 
and  third  order  mixing  properties.  In  addition,  incommensurate 
chains  have  been  seen  to  exhibit  complete  mode  locking  over  the 
entire  range  of  dc  fields  and  external  frequencies.  (These 
results  resolve  a  dispute  between  the  Illinois  and  the  AT&T  Bell 
Laboratories  groups.  The  manuscript  is  being  prepared  for 
Physical  Review  Letters  and  the  abstract  is  included  in  section 


Ill .  Publications 


1.  Sliding  Dynamics  of  the  Incommensurate  Chain 
Physical  Review  Letters,  52.,  65  (1984). 

2.  Dynamics  of  Incommensurate  Structures:  An  Exact 
Solution 

Physical  Review  (Rapid  Communications)  ,  B30 .  2974 

(1984)  . 

3.  Electromechanical  Properties  of  Charge  Density  Wave 
Conductors 

Physical  Review  Letters,  56,  1194  (1986) . 

4.  Oscillatory  Instability  in  the  Dynamics  of 
Incommensurate  Structures 

Physical  Review  Letters,  58,  1903  (1987) . 

5.  Mixing  in  Charge  Density  Wave  Conductors 

Physical  Review  (Rapid  Communications)  B35 ,  7745 

(1987)  . 

6.  Mode  Locking  in  Charge  Density  Waves:  a  Classical 
Theory 

to  be  submitted  to  Physical  Review  Letters. 


v.mi  Mh  ij.  N.  mhik  I  PHYSICAL  REVIEW  L  L 1  1  I.  R  >  2  J»m  am  l'W( 


Sliding  Dynamics  of  the  Incommensurate  Chain 

Leigh  Sneddon 

Martin  Fisi’fr  School  or  Ph\ sics ,  Brandcis  l'niversit\ .  Waltham.  Massachusetts  01125! 

(Received  S  July  19*3) 

The  dc  dynamics  of  the  sliding  incommensurate  chain  is  reduced  to  a  purely  static  prob¬ 
lem.  The  sliding  system  is  described  by  a  static  hull  function  which  becomes  singular, 
above  the  critical  pinning  strength,  as  the  velocity  approaches  zero.  Both  ac  and  dc 
sliding  dynamics  are  determined  numerically  for  the  cases  of  weak  and  strong  pinning 
and  short-  and  long-range  interactions.  Excellent  agreement  is  obtained  with  experiments 
on  sliding  charge-density  waves  near  threshold,  both  in  N’bSej  and  in  TaSj. 

P ACS  numbers;  72.15.Nj,  72.15.Eb,  72.20. Ht 


Incommensurate  structures,  such  as  charge- 
density  waves  (CDW’s)  and  adsorbed  monolayers, 
are  now  familiar  in  solid-state  physics.1  The  dis- 
covery:  of  electrical  conduction  due  to  sliding 
CDW's  raised  a  wide  range  of  questions  concern¬ 
ing  the  dynamics  of  sliding  incommensurate  struc¬ 
tures.  This  article  reports  progress  in  the  ana¬ 
lytic  study  of  such  dynamics  and  in  the  under¬ 
standing  of  related  experimental  results. 

The  incommensurate  system  studied  here  is  a 
simple  extension  of  the  model  of  Frenkel  and 
Kontorova,1  and  the  dimensionless  equations  of 
motion  can  be  written 

O  j  Yj  p  F  p  (2  U  j  ~  0  j  -  p  ~  tt  j  +  p  ) 

-  f  *■  PcosH  (j  +  Hj ),  (1) 

where  j=  1,2 . .V,  P  is  the  strength  of  the  pin¬ 

ning  force,  and  2tr ///  its  wavelength.  The  case  oL 

v-f-PF.lu  J  , 

w^-(2\l  ~Yp  A’pcosp/y^l  -  ipJlv)w  ^  *■  PF^wJ 
for  1  ti  '  X'2,  where 

V 

FA't’Ji  £.V  Tj  <i  cos tl\j  -  2  Y  n.cosl 

1  - 1  1  -  m  -  N’  1 

'Hileficwirnl  Dt  tvv 

For  constant  v,  static  solutions,  with  /V^-0 
for  1,  were  found.  A  finite  number. 
of  the  m  _  were  retained  and  the  corresponding 
coupled  equations  Mb)  were  solved  numencallv, 
treating  the  retained  a  t  to  all  orders.  The  coup¬ 
lings  in  Me  i  were  computed  with  V  large  enough 
that  increasing  it  had  no  significant  effect.  These 
solutions  were  tested  m  two  wavs.  Firstly,  the 
dp  characteristic  was  calculated.  With  weak  pin¬ 
ning.  clear  convergence,  with  increasing  , 
to  a  linear  response  at  t  0  wa.s  found,  with  :v> 
thre-'r.ol  I.  For  strong  pinning  'lie  emergence  ->t 
(  •hr>'-,hnld  singularity,  with  increasing 


X—x>  and  ///2tr  equal  to  an  irrational,  p,  was 
studied  by  considering  z//2g=  Af  v '.V,  where  M s 
is  an  integer,  having  no  common  factor  with  .V, 
and  Yfv/,V— p  as  X—  «. 

Equation  (1)  was  Fourier  transformed  to  re¬ 
place  the  w/s  by  phase-shifted  Fourier  compo¬ 
nents  wm  =  a„  expUoJ: 

N 

wm  =exp(- imtfw^X'1  Y  exp(- imHj)uj  (2) 

/=i 

for  m  -0, 1,  2, ... ,  .V-  1.  Using  Bessel  functions 
to  expand  the  cosine  in  (1)  gives  terms  with  an 
explicit  X  dependence.  As  .V—®  for  a  bulk  veloc¬ 
ity  v  >  0,  these  terms  can  be  shown  to  vanish 
to  every  finite  order  of  perturbation  in  P  and 
oscillate  at  arbitrarily  high  frequencies  about 
zero.  To  treat  the  limit  X—  »,  only  those  terms 
with  no  explicit  ,V  dependence  were  retained. 

This  gives 

(3a) 

(3b) 

/U  (3c) 


was  clearly  indicated.  The  results  thus  agree 
with  present  knowledge1  at  low  velocities  and  are 
also  correct  to  all  orders  in  perturbation  theory 
at  moderate  and  high  velocities.  Secondly,  the 
solutions  were  tested  for  stability  to  small  per¬ 
turbations  and  were  found  to  be  dynamically 
stable. 

Since  these  solutions  are  static,  exploiting 
translational  invariance  has  transformed  the  dc 
sliding  l\  n.i  tines  d  the  infinite  incommensurate 
chain  '  o  a  pur<-!v  static  problem. 

From  2\  th'1  r>u ric r  component s.  \  of  the 
'list  .rtiotis  Mien  have  a  simple  lime  dependence 


h> 


'  F'S  l  [  he  Viper  i. an  lVrok.il  Sm; 
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>•„  =  const  x  corresponding  result  for 

excitations  about  the  v  =  0  limit  was  obtained  by 
Novaco.M  Further,  there  is  a  static  hull  function 
describing  dc  sliding  solutions  with  v  >  0.  That 
is,  u  j  -  vt  is  a  function  of  a  single  variable: 

u,  -  it  =gt(j  -  vt),  (4) 

where  the  periodic  function  g  is  given  by 

gv{x)  =  2  2  a„cos(m//x  *  a  J.  (5) 

l  -  m<  N  J2 

Much  attention  has  been  given4,6  to  the  way,  at 
f  -  v  -  0,  gv  changes  from  being  analytic,  for  weak 
pinning,  to  singular,  for  strong  pinning.  The 
present  studies  showed  the  amplitudes  am  de¬ 
caying  exponentially  with  m  when  r>0,  even  with 
strong  pinning.  The  function  gv  is  then  analytic 
for  t  >  0.  Thus  for  strong  pinning  a  breaking-of- 
analyticity  transition  occurs  in  the  new  hull  func¬ 
tion,  as  the  velocity  approaches  zero  at  thresh¬ 
old.  The  complicated  time  dependence  of  the 
Uj(t)  near  threshold  is  expressed,  by  (4),  com¬ 
pletely  in  terms  of  the  emergence  of  singulari¬ 
ties  in  this  new  hull  function. 

Linear  ac  response,  in  the  presence  of  a  dc 
field,  has  been  studied  experimentally7  in  the 
CDW  systems  NbSe,  and  TaS,.  At  fields  a  few 
times  threshold  low-order  perturbation  theory3 
is  not  useful;  but  this  region  is  the  most  common¬ 
ly  studied  experimentally  because  the  nonlinear 
effects  are  larger  than  in  the  high-field  region, 
and  sample  heating  is  not  a  problem.  The  pres- 


0  2  04  06 
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ent  techniques  were  therefore  used  to  determine 
the  ac  response  of  the  sliding  incommensurate 
chain  near  threshold.  Having  reduced  the  dc 
dynamics  to  a  static  problem  is  very'  useful.  The 
ac  response  of  a  time-dependent  solution  is  much 
more  difficult  to  obtain  than  that  of  a  static  solu¬ 
tion.  for  which  the  ac  response,  like  the  dc  char¬ 
acteristics,  can  be  determined  without  numerical 
integration. 

The  CDW’s  in  NbSe,  and  TaS3  are  three-dimen- 
Sionally  coherent.'’  One  effect  of  higher  dimen¬ 
sionality  is  to  increase  the  coordination  of  the 
system.  To  mimic  this  increased  coordination 
crudely,  a  sixfold-coordinated  chain  was  con¬ 
sidered  With  =  A'=  =  A,  =-j-;  Kp  =  0,  p  >  3. 

By  considering  in  (3)  a  small  perturbation  about 
a  static  dc  solution,  the  ac  response,  cr(uj)  =  o' 

+  its"  was  determined,  for‘°^)=  (51/2  +■  l)/2  and 
HP- 3.0.  The  results  (with  4mlI  =  15)  for  a'  and 
the  dielectric  response  -  o"/u>  are  shown  in  Fig. 
1.  The  basic  features  in  Figs.  1,  3,  and  5  are 
preserved  with  increasing  Umax-  The  threshold 
force  was  estimated  from  the  dc  results. 

Figure  2  shows  experimental  results7  for 
Recta)  and  eta)  of  the  sliding  charge-density 
wave  in  NbSe3.  Figure  1  is  seen  to  account  well 
for  the  voltage  and  frequency  dependence  of  both 
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FIG.  3.  ac  response  of  incommensurate  chain  with 
infinite- range  interactions;  c.f.  Fig.  4  and  text. 


components  of  the  ac  response.  This  may  not 
have  been  expected  since  CDW  dynamics  are 
dominated  by  randomly  positioned  defects11  while 
the  chain  is  in  a  periodic  potential. 

In  experiments7  performed  on  TaS3  at  130  K, 
the  sharp  interference  features  seen7  with  NbSe3 
(Figs.  1  and  2)  were  not  observed.  TaS3  becomes 
a  semiconductor12  below  the  CDW  transition, 
while  XbSe3  is  metallic.2  At  130  K  the  conductivi¬ 
ty  of  TaS3  has  fallen  two  orders  of  magnitude 
from  its  value  at  tne  transition.  As  discussed 
earlier, 13  this  reduces  the  screening  capacity 
of  the  normal  electrons  and  can  allow  long-range 
Coulomb  interactions  of  the  CDW  with  itself. 

The  sliding  dynamics  of  Eq.  (1)  with  long-range 
interactions,  Kp  =2/.V  for  all  p,  was  therefore 


determined.  The  results  (with  p  . 


=  20)  are 


shown  in  Fig.  3,  and  can  be  compared  with  the 
experimental  results  in  Fig.  4.  Not  only  does 
including  long-range  interactions  account  for  the 
absence  of  interference  features,  but  the  proper¬ 
ties  of  the  incommensurate  chain  are  seen  to 
match  those  of  TaS3  extremely  well.  The  dif¬ 
ference  between  the  ac  properties  of  NbSe.,  and 
TaS3  can  now  lie  understood  for  the  first  time  as 
being  due  to  the  presence  in  TaS3,  as  suggested 
earlier,13  of  long-range  Coulomb  interactions 
of  the  CDW  with  itself. 

The  ac  response  was  also  determined  with  r 
=  0,  and  compared  to  the  dc  conductivity  '  t . 

The  results  'with  p..  =20)  are  shown  in  Fig.  5 

for  long-range  interactions.  Similar  results 
were  obtained  for  the  sixfold-coordinated  chain. 
The  experimentally  observed"  •  15  scaling,  of 
field-  and  frequency-dependent  conductivities. 


bias  voltage  Vdc  (mV) 

FIG.  4.  ac  response  (fief.  7)  of  TaS3. 


is  thus  exhibited  by  this  classical  model,  and  can 
no  longer  be  regarded15  as  evidence  for  a  quantum 
mechanical  theory  of  CDW  conductivity. 

Although  these  calculations  do  not  probe  asymp¬ 
totic  low-frequency  threshold  properties  it  is 
interesting  to  speculate  that  the  detailed  form  of 
the  potential  becomes  less  important  as  one  ap¬ 
proaches  threshold.  In  any  case,  the  comparison 
of  theory’  with  experiment  seen  in  Figs.  1-5  shows 
that,  in  fields  comparable  to  threshold,  the  in¬ 
commensurate  chain  gives  a  much  better  pic¬ 
ture  of  CDW  dynamics  than  might  have  been  sus¬ 
pected. 
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FIG.  f),  Scaling  of  licld-dcpcndont  (crosses)  and 
freqnenev-dependent  (circles)  conductivities  in  (he 
classical  chain. 
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A  simple  exact  solution  of  the  incommensurate  chain  with  infinite-range  interactions  is  obtained  with  the 
use  of  both  analytic  and  graphical  techniques.  The  ground  slates  and  all  metastable  states  are  identified. 
Results  are  derived  for  the  deptnning  transition,  the  sliding  threshold,  and  the  excitation  spectrum  as  a 
function  of  pinning  strength  and  applied  field  The  ac  response  is  also  obtained  and  has  a  low-frequency 
singularity  at  threshold,  but  the  dielectric  constant  is  bounded,  as  seen  in  charge-densitv-wave  experiments 


Incommensurate  structures,  such  as  charge-density  waves 
(CDW's)  and  adsorbed  monolayers,  are  well  known  in 
solid-state  physics.1  The  incommensurate  harmonic  chain 
has  recently  been  seen2  to  provide  a  useful  model  of  CDW 
conductivity  and  long-range  interactions  have  been  seen2,3 
to  be  important  in  semiconductor  CDW  systems,  e.g.,  TaSj. 

If  every  particle  in  a  chain  interacts  equally  with  every 
other,  the  equation  of  motion  can  be  written 

— =  j  ( t)-r  P  sin  (  Hi  +  Uj)+  (  U)  —  Uj  .  (1) 

Here  j  —  i ,  2 . 

,\ 

<  U)  =  lim  .V" 1  V  u  , 
v  — 

P  is  the  strength  of  the  pinning  force,  and  2ni  H  its 
wavelength,  an  irrational.  Fisher4  has  studied  the  related 
problem  where  P  is  replaced  by  a  randomly  distributed  vari¬ 
able.  The  purposes  of  this  Rapid  Communication  are  to 
provide  a  considerably  simpler  exact  solution  of  the  incom¬ 
mensurate  case  (fixed  P)  and  to  derive  a  number  of  new 
exact  results. 

Solutions  to  U)  can  be  written 

Lij  ( II  —  n  I  1 1  —  g  (  Uj  —  (i  ( t ) ;  t )  ( 2 ) 

where1  a(  x  —  2  — ; / 1  =  ?  (  x.t ) .  and 

/  ( t  i  =  -  I  2  -  )  ~  1  J  dx  y  (  x , : )  ( 3 ) 

Substituting  1 2)  in  f  I )  and  using  the  identity 

1 1 m  V  ~ !  'V  /•  ( y  )  )  =  r  p(  x  >  tlx 

TT,  -■ 

for  y  irrational  and  I  I  x  —  1  )  =  / '(  x  )  gives 


rill 


-  [VL 

hx 


=  P  sint  x  -  y )  -  g  . 
or 


where,  for  bulk  velocity  d/l-iii  L'i,Jt.  satisfies 

i  d  1 1  -  Mr) 


( 4 ) 


15) 


Thus,  ’he  form  *2)  yields  an  equation  of  motion  which 
depends  an  >  a ,/  <  t'i/Jt.  hut  not  on  <  U > .  exhibiting  the 
■ransi.uion.il  invariance  ot  the  incommensurate  system.'’ 

I  he  T.insient  and  hvstcrctic  properties  of  (|)  are  current- 
v  inner  .nvestat.iuon  The  rest  ot  this  article  will  consider 


constant  /,  in  which  case  a  and  v  are  constant,  and  one  can 
search  for  solutions7  with  0g/9r  =  O.  Equations  (2).  (4), 
and  (5)  then  simplify8  to 

Uj(  i )  =a  -t-  g{  Hj  -b  a  )  ,  1 6) 

v(  1  +  g')  =  P  sin(x  +  g  )  —  g  .  t7) 

When  u  =  0.  (6)  is  a  simple  transcendental  equation  which 
can  be  solved  graphically.  For  P  <  1  here  is  a  unique  solu¬ 
tion  [Fig.  1(a)].  It  is  continuous  and  odd  so  that,  using  (3), 
when  u  =  0.  /  =  0,  and  there  is  no  sliding  threshold.  For 
P  >  1  there  are  multiple,  discontinuous  solutions,  g,  many 
with  nonzero  means  [Figs.  Kb)  and  1(c)],  The  threshold 
force  is  clearly 

fT=  max  —  ( 27t )  ~ 1  f  g{x)dx  .  (8) 

111  J 

Thus  the  critical  value  of  P  defining  the  deptnning  transition 


FIG  I  Stationary  -talcs  1  j  •  —  '  b '  /’ -  L  ground 

states  lor  '  —  0.  h),  n  <  ■  -  ,  ' 11 1  ■  '  ”  ' '  ■'  *■  ty.  in)  <ci  P  -  \ 

a  meiasiable  slate 
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below  which  the  sliding  threshold  and  multiplicity  of  solu¬ 
tions  disappear,9  is  immediately  seen  to  be  Pc  =  1. 

The  threshold  force  can  be  determined  graphically,  using 
(8),  for  all  P.  Further,  in  the  limit  P —  Pc+ ,  the  mul¬ 
tivalued  region  in  Fig.  1  shrinks  towards  the  origin  where 
(7)  gives 

x-g3/6+(P- Pc)g  =  0  ,  (9) 

and  (8)  gives,  with  elementary  calculus,  J\=  ( 9/ < 4rr ) ] (P 
-  Pc)*t  as  P  —  Pc+,  where  Tr=  2.  For  large  P  one  obtains 
fr- P~TT  +  0(P~yn). 

Further,  it  is  immediately  clear  from  Fig.  1(b)  that  at 
/—  Jt  there  is  only  one  stationary  state,4  gT- 

Turning  to  dc  dynamics,  a  =  vt  in  (6)  and  g  is  continuous. 
As  v—  0,  g  will  approach  gT  as  / —  ff  ■  It  is  easy  to  see 
from  (7)  that  for  x  away  from  the  critical  value  xT , 
g  -  gr=0(v).  In  the  vicinity  of  xT ,  putting  x  =  xT+y, 
x  4-  g(x)  —  xr+  gT^xf  )  +  h  (>■ )  and  considering  (7)  in  the 
limit  of  small  y  and  h(y )  gives10  vdh/ dy  =  y  +  ah2,  where 
a  =  —  gr(xf  )/2.  Transforming  by 

—  ah(/3z)w(z)  ,  (10) 

dz 

where  a  —  (a2/v)l,}  and  /3 ■=  (v2/ a  )l/}  gives  w"=  —  zw(z), 
the  solutions  of  which11  are  the  Airy  functions  Ai(  -  r)  and 
Bi(  —  r).  The  boundary  condition  that  h(y)  be  finite  and 
negative  for  y  <  0  as  v—  0  means,  integrating  (10),  that 
w(z)~  0  as  c  —  —  oo.  This  eliminates  Bi(  —  z),  which 
diverges  as  r—  Thus 

l/  1/3  Ai'(  -  (a/v2)U3y) 

a 2  Ai(  —  (a/v2)l,3y ) 

and  the  limiting  value  of  y  =  x-xr  for  finite  h,  as  v—  0,  is 
( v2/a)inzo  where  z0  is  the  first  zero  of  Ai(-z).  This 
result  is  seen  graphically  to  give  a  dominant  contribution 
~v2/3  to  f—fr  in  (3).  Thus  v  =  B(f - /r)3/2,  where 

[27ra1/3/(zoA)l3/2  and  A  =  gT(xf  )  —  gr(xf )  ■  Further, 
as  P  —  P*  using  (9)  gives 

B  =  ( 2,/3rr/ ( 3z0)  13/2(P  —  Pc) _l/2  . 

Thus  the  depinning  transition,  and  the  y  threshold  ex¬ 
ponent  with  a  coefficient  which  diverges  as  P  —  Pc+,  in 
agreement  with  Ref.  4,  can  be  obtained  quite  straightfor¬ 
wardly  for  the  incommensurate  chain.  The  simplicity  of  the 
present  calculation  is  due  to  the  elimination  of  the  center- 
of-mass  coordinate  (  U)  leaving  a  function  of  only  one  vari¬ 
able  to  be  determined  in  solving  incommensurate  dc  dynam¬ 
ics. 

Further,  it  is  possible  to  determine  the  energy  and  stabili¬ 
ty  of  each  stationary  state  and  thus  specify,  for  P  >  Pf, 
which  is  the  ground  state,  which  are  the  metastable  states, 
and  which  are  the  unstable  states.  The  energy  correspond¬ 
ing  to  Eq.  (1)  is 

H-^PcosiHj  +Uj)  +  (4N)-'X<Ui-Uj)2  ■ 
i  u 

Using  16)  and  choosing  g  to  minimize  H  shows  that  the 
ground  state  has  a  single  discontinuity,  which  moves  from  0 
to  xT  as  /  increases  from  0  to  JT  (Fig.  1(b)). 

To  distinguish  metastable  from  unstable  states  the  linear 
stability  of  stationary  solutions  \U/\  of  (1)  must  be  studied. 
Replacing  Ut  by  Ut  +  itj  in  (1)  gives,  to  OlUj).u 


=  +  Xf\)u,  where  \t&u>  =  (PcosI  Hj  ^  U, )  -  1  )o7  anu 

.Vflu>  =  N~[.  The  Green's  function  £K  r  )  =  ( :  -  A/0 
-/Wi)-1  is  easily  calculable  by  summing  the  usual  expan¬ 
sion  in  powers  of  Tf,  With 

Gou>  (z)  =  ( z  —  Kj)~'hfj,  \j=  Pcosi  Hj  U,)  —  1  . 

one  obtains  G  =  Gnf  GnTGn.  where  TiJ)  =  ( ,V(  1  -  /a  )  ]  “ 1 
and  |iz(c)  =  A_1V/(z-A/)_l.  Thus  the  poles  of  G,  which 
are  eigenvalues  of  the  stability  problem,  are  the  roots  of 
/Kz)=  1.  A  plot  of  *K\)  for  real  \  then  immediately 
shows10  that  any  g  which  occupies  a  finite  part  of  the  middle 
branch  (dashed  line  in  Fig.  1)  will  be  unstable.  Further,  as 
/V—  oo.  the  largest  root  of  m(-)“1  remains  an  isolated  ex¬ 
citation.  Transforming  sums  to  integrals,  it  is  straightfor¬ 
ward  to  show  that  this  eigenvalue  is  negative  for  any  state  g 
which  avoids  the  “unstable”  middle  branch.  Thus  any  such 
state  which  is  not  the  ground  state  is  metastable.  An  exam¬ 
ple  of  a  metastable  state  is  shown  in  Fig.  1  (c). 

Since  the  Green’s  function  is  known,  the  exact  density  of 
relaxationai  1(1)  is  massless]  excitations.  p<\)  can  be  deter¬ 
mined  in  the  usual  way.  For  P  <  Pc,  transforming  sums  to 
integrals,  and  writing  p(X)  =  - x(P2  -  (X  - <- 1  ):)-i/2/t\  one 
obtains  for  large  N 

p(X)  =  p(X)0(P2— (X+  IHJ  +  .V-'S'X)  , 

which  is  shown  in  Fig.  2(a).  The  isolated  excitation  at  x  =  0 
is  the  sliding  mode  of  the  unpinned  chain. 

For  P  >  Pc,  /  =  0,  the  excitation  spectrum  p(x)  was 
determined  for  the  ground  state  shown  as  curve  (i)  in  Fig. 
1  (b).  The  result  for  large  .V  is 

p(X)«p(X)0(l  +  P  +  X)0(Xo-X)  +  iV-’SU-M  ■ 

where  X0=  -  1  +  Pcoslg(0+ ))  <  0,  and  x0  <  X/>  <  0.  This 
spectrum  is  shown  in  Fig.  2(b).  As  P—  P*.  (9)  gives 
X0=  -  2(P  —  Pc).  Further,  as  P—P*.  XP=-MP-PC) 
where  b  is  between  0  and  2  and  is  the  root  of 


-2  -I  0 


-4  -3  -2  -4  -3 


FIG.  2.  Excitation  spectra  fi(x)  (al  /’-I)  5.  / -0,  ' h )  /’ -3. 
/  -  0;  (c)  P  —  3,  0  <  /  <  fT. 
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Thus  as  P  —  P*  the  gap  in  the  ground-state  excitation 
spectrum  doses  as  iP-  Pc)  ’  where  ii|  =  1. 

For  P  >  Pc ,  the  excitation  spectrum  was  also  determined 
for  f>  0.  For  the  ground  state  shown  as  curve  (ii)  in  Fig. 
I (b),  one  finds 


p(A).  -  1  -  P  <  A  <  A*  . 


p  (  A  )  = 


p(  A  )/2.  A  *  <  A  <  A  _  . 
V-  'S (  A  —  A  f ) .  A_  <  A  . 


where  A  *  =  -  1  -r  P  cosl  .xy  g  (  .xy1 )  1 .  and  A  _  <  \f  <  0.  As 
/-  ff,  A  _(/')= -.9  (yr-/)1/:.  where  .4  ( /> )  —  ( 2/3  )1/2 
as  P —  /V  This  spectrum  is  shown  in  Fig.  2(c).  Further, 
as 


/—ff.  k,  =  \  -  ~  D  exp(  —  C ( J'T  —  f)~U2]  , 


FIG.  3.  dc  characteristic  for  P—  1x5  (o):  also  (wuh  matched 
thresholds  and  targe  v  slopes)  for  single-particle  model  (Ref.  12) 
(+  )  and  NbSe3  CDW  (Ref.  16)  (A). 


where  D  -  2(  P2  -  i )  P  and  C  =  ( P2  -  1  )/(4a  )  ]1/2.  That 
is.  as  /—  ff ,  the  isolated  excitation  collapses  exponentially 
rapidly  onto  the  edge  of  the  continuum.  Moreover,  the  gap 
in  the  excitation  spectrum  closes  as  (/r  —  f)*1  where 

'it  2  "  y . 

Finally,  knowing  the  Green’s  function  gives  the  exact 

- I 


linear  ac  response  of  the  pinned  lattice.  For  P  >  Pc  the 
result  is 

rr<<o,f  \P)  =  ito  jl  —  2ir(A/  —  i<a(  E(  loi.  A_  ) 

+  £(/cu.A+))]-'|  .  (11) 

where  A f*-g(xf)-g(x/~ )  and 


E  (:,Af 


2 


[P2-(\-:)2Yn 


coth" 


P+\~: 


P  —  1  —  A 


[p2_(l_.)2]l/2  (/.+  1  +  X]l/2 


The  result  (11)  may  be  compared  with  CDW  experiments 
for  both  the  real  and  imaginary  parts  of  or  (a >)  for  all 
i/I  *£  Jr  and  for  all  <u.  The  single-particle  model15  is  easily 
seen  to  have  the  property  lim._0- Imo7oj  —  oo  as  / '  —  fT. 
No  such  divergence  has  been  observed  experimentally.13  As 
id  —  0  and  f  —  fr.  the  term  in  square  brackets  in  (11)  be¬ 
comes 

P(|A_I  —  mq) 

2  ( P5  -  1 ) 

While  !a_  |  —  0  as  / —  />,  the  dielectric  constant,  t  —  kt/ui. 
nevertheless  remains  finite  near  fT  even  as  cu  — ■  0.  The  en¬ 
tire  spectral  weight  of  the  single-particle  model  is  at  one  fre¬ 
quency  which  approaches  zero  at  threshold.  The  pinned 
many-body  system  considered  here  has  a  broad  spectrum  of 
excitations,  only  one  edge  of  which  approaches  zero  at 
threshold.  This  is  the  first  explicit  demonstration  that  in 
such  a  case  the  singularity  in  n  is  weaker  than  in  the 


I 

single-particle  case  and  e  can  remain  finite,  as  is  observed 
experimentally.13 

For  nearest-neighbor  interactions,  the  ground  state  g  is 
more  singular  and  some  numerical  determinations  of  depin- 
ning  exponents  for  a  stationary  center  of  mass  have  been 
made.9  14  Equation  (1)  can  be  considered  a  mean-field 
theory  for  the  incommensurate  chain  with  finite-range  in¬ 
teractions.  The  present  straightforward  solution  may  be 
useful  in  the  ultimate  analytic  solution  of  sliding  threshold 
dynamics  with  short-range  interactions.  For  constant  /  >  fT 
outside  the  threshold  region,  the  truncation  procedure  of 
Ref.  2  gives  an  accurate  solution  of  Eq.  (1),  including  the 
response  to  small  ac  perturbations.  As  seen  for  example15 
in  Fig.  3,  this  combines  with  the  present  results  for  /  <  fT 
to  give  a  complete  steady-state  solution  of  this  system. 

The  author  thanks  E.  Gross  for  stimulating  discussions 
and  the  Air  Force  Office  of  Scientific  Research  for  support 
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A  model  of  two  mutally  incommensurate,  interacting,  dynamical,  many-body  systems  is  present¬ 
ed  and  solved.  Quasiperiodic  forms  are  shown  to  describe  both  the  dc  properties  and  the  complete 
set  of  linear  excitations.  By  use  of  one  system  to  represent  the  crystal  lattice  and  the  other  a 
charge-density  wave  (CDW),  all  the  recently  discovered,  and  as  yet  unexplained,  electromechanical 
properties  of  CDW  conductors  are  shown  to  occur  in  this  model.  The  internal  degrees  of  freedom 
of  the  CDW  are  shown  to  be  of  central  importance. 

PACS  numbers:  72.15.Nj,  72.15.Eb,  72.20.Ht 


Until  recently  the  widespread  interest  in  charge- 
density-wave  (CDW)  conductors  has  centered  on  their 
nonlinear  electrical  properties.1  Elegant  experiments 
by  Brill  and  Roark2  and  by  Mozurkewich  et  al.?  how¬ 
ever,  have  now  shown  that  the  motion  of  a  CDW  also 
changes  the  mechanical  properties  of  the  host  crystal. 
Speculation  as  to  the  origin  of  these  phenomena  has 
included  CDW  domain-wall  relaxation2  and  phason- 
phonon  interactions3;  and  a  satisfactory  calculation  of 
the  current  dependence  of  mechanical  properties  is 
still  lacking. 

This  Letter  addresses  the  question  of  whether 
single-domain,  noncommensurate,  classical,  many- 
body  models,  which  have  exhibited  a  variety  of  non¬ 


linear  electrical  properties  of  CDW’sw  (e.g.,  dc 
characteristics,  ac  response,  ac-dc  interference,  and 
voltage  fluctuations)  can  also  give  a  satisfactory  ac¬ 
count  of  these  new  electromechanical  properties.  A 
model  is  solved  which  includes  both  the  host-crystal 
degrees  of  freedom  and  the  internal  degrees  of  free¬ 
dom  (IDF’s)  of  the  CDW.  It  consists  of  two  interact¬ 
ing  mutually  incommensurate  chains,  one  representing 
the  crystal  lattice  and  the  other  the  CDW.  Conserva¬ 
tive  and  viscous  interactions  between  the  chains  are 
used  to  represent  respectively  the  pinning  force  which 
produces  the  CDW  threshold,  and  the  dissipation 
which  results  from  CDW  motion  through  the  crystal. 

Dimensionless  classical  equations  of  motion  for  this 
model  can  be  written: 


"1* a  +  2*Wa-e+P-Xj  mxat)(*a -  Uj)  =  ~n 2jF(Xaj)  -  El.  (1) 

2,D,  Uj  +  la  ^  Xaj  >  < (  Uj  -  <t>a  )  =  2a  F( XaJ  )  +  E.  (2) 


where  m  is  the  ionic  mass,  the  inertia  of  the  CDW  is 
negligible.  (t>a  and  Uj  are  the  displacements  of  particle 
a  in  the  crystal  lattice  and  particle  j  in  the  CDW, 
respectively,  and  2na  and  Hj  are  their  respective 
undisplaced  positions,  so  that 

Xa]^h Ta  +  *a-[Hj  +  Uj)  (3) 

is  the  distance  between  them.  The  coefficients  A3  and 
Dp  are  the  spring  constants  of  internal,  harmonic  re¬ 
storing  forces  in  the  lattice  and  the  CDW,  respectively. 
The  function  W(x)  is  a  weighting  function  centered 
at  .V  “  0,  and  represents  the  spatial  range  of  the  dissipa¬ 
tive  interactions.  The  force  F between  particles  a  and 
j  also  depends  on  their  separation  Xaj.  E is  the  electric 
field  acting  on  the  CDW  and  -  EL  is  the  force  which 
keeps  the  lattice  stationary.  The  incommensurate  limit 
is  approached  by  a  consideration  of  M  particles  in  the 
crystal  lattice  chain.  .V  in  the  CDW  chain,  and  Ml  N 
=*  HUtt  —  a  fixed  irrational.  The  parameter  /x  is  a 
formal  expansion  parameter  ultimately  to  be  set  equal 
to  1.  This  model  thus  treats  the  lattice  and  the  CDW 
as  two  separate  entities,  each  distorting  the  other.7  In 


the  stiff-lattice  limit  (A  »  F),  <t>a=  (<t>ct).8  and  the 
crystal  lattice  simply  produces  a  rigid  incommensurate 
potential  through  which  the  CDW  moves,  exhibiting  a 
variety  of  experimentally  observed  nonlinear  electrical 
properties.4*6 

The  first  general  result,  which  can  be  checked  by 
substitution  for  any  A,  is  that  stationary  states,  and 
states  of  dc  relative  motion  of  the  centers  of  mass  of 
the  two  interacting  chains,  are  described  by  quasi¬ 
periodic  forms 

=  Q(y).  Uj(t)  =  vt  +  G(x), 

y  =  2na  —  vt.  x^Hj  +  vt,  (4) 

Q(y  +  H)-Q(y).  G(  x  +  2tt  )  -  G(x ), 

where  v=(Uj),  the  CDW  center-of-mass  velocity. 
Thus  the  quasiperiodic  form  first  shown  to  describe 
the  statics  of  a  single  incommensurate  chain.9  and  sub¬ 
sequently  shown  to  describe  its  dc  dynamics  as  well,4  is 
now  seen  to  describe  both  the  statics  and  dc  dynamics 
of  two  interacting  dynamical  systems. 
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Further,  the  linear  fluctuations  of  the  dual  interact¬ 
ing  system,  both  when  static  and  when  undergoing  dc 
relative  motion,  are  likewise  characterized  by  quasi- 
periodic  forms.  If  <ba  =  <ba  exp(  -  iut )  and  Uj  =  Uj 
*exp(  —  /co/ )  are  fluctuations  in  the  crystal  lattice  and 
CDW  chains,  respectively,  a  complete  set  of  linear 
fluctuations  is  given  by 

<ba(  i )  =  e“,2war(  y ).  Uj(  i)  =  e^hix), 

(5) 

r (>■  +  H)  =>  /(>'),  h  ( r  +  2tt  )  =  h  ( x ) . 

where  a  and  y  are  as  in  (4).  The  normal  modes  of  the 
dual  interacting  system  thus  resemble  Bloch  waves.10 
If  the  two  chains  do  not  interact,  r  and  h  are  constant 
and  the  excitations  are  simple  traveling  waves.  Inter¬ 
chain  interactions  produce  a  quasiperiodic  modulation 
of  each  wave.  In  the  special  case  of  the  stiff-lattice 
limit  these  excitations  reduce  to  those  of  incommensu¬ 
rate  single-chain  systems,  which  have  been  studied  nu¬ 
merically  in  the  static  case  (u  =  0)9:  We  now  see  that 
these  excitations  are  “quasiperiodic  Bloch  waves.” 
Further,  the  exponential  factors  in  <t>a  and  Uj  are  waves 
with  phase  velocities,  respectively,  of  io/q  in  the  crystal 
lattice  and  <n/q  —  v  in  the  CDW.  Since  the  CDW  is 
traveling  at  v  relative  to  the  crystal,  both  waves  are 
traveling  at  the  same  velocity  in  the  laboratory  frame. 

To  compare  the  properties  of  this  model  with  exper¬ 
imental  data,  a  small  force  will  be  considered 

to  act  on  the  crystal,  and  the  resonant  frequency,  a>o, 
defined  as  the  frequency  at  which  the  driving  and  the 
displacement  are  tt/2  out  of  phase,2  will  be  determined 
as  a  function  of  the  voltage  across  the  system.  The  in¬ 


verse  displacement  amplitude  at  resonance  is  then  a 
measure2  of  “internal  friction,”  Q~l 

Since  the  effect  of  the  CDW  on  the  Young's 
modulus  is  of  the  order  of  a  percent,2  Eqs.  (1)  and  (2) 
will  be  studied  as  an  expansion  about  the  stiff-lattice 
limit.  This  is  most  clearly  done  by  expansion  of  the 
solution  in  powers  of  the  formal  parameter  m.  which 
premultiplies  the  forces  distorting  the  lattice. 

The  calculation  begins  by_  replacement  of  Uj  by 
Uj  +  Uj(t)e~ 4>q  by  <Pa+<ba(t)e~ and  -  EL  by 
—  El  +  fae~iQ“  in  Eqs.  (1)  and  (2),  of  which  the  U/s 
and  <t>a’s  form  a  dc  solution;  retention  of  only  terms 
linear  in  w,  <t>,  and  /;  and  consideration  of  6a  of  the 
form,  to  Oifx °),  of  a  normal  mode,  i.e., 

4>a  =  exp(  iqlna  ).  (6) 

Multiplied  throughout  by  exp(  -  iq2na),  and  aver¬ 
aged  over  a,  Eq.  (1)  then  gives 

/+  mw2  —  A ( // )  =  bf  (7) 

where11 

ZaJ  =  [F(XaJ)-i«>W(Xaj)}(\-Uje-*2"a).  (8) 

To  determine  /  to  0(/x)  in  (7)  only  requires  that 
Z„y,  and  hence  f/;,  <t>„,  and  uJy  be  determined  to 
o(n°),  i.e.,  in  the  rigid-crystal-lattice  limit,  for  which 
we  may  put  d>a=  (<J>a)  =0.  To  determine  Uj  and 
the  case  of  infinite-range  internal  CDW  interactions5 
(“mean  field  theorv”)  will  be  considered: 
ZpDpUj  p~  Uj-(U).  Use  of  Eqs.  (3)-(6)  and 
linearizing  of  (2)  then  gives  (for  q^O  so  that 

( uj)  0) 


v  W<  Xa) h'( x )  +  [  1  +  Xu !  C,  +  iqv  W(  Xa) }  ] *  (x )  =  2<»Ca  exp[ iq ( 2ira  - x)  ],  (9) 

where  Xa  =  2 ira  —  x  —  g(x),  C„  =  F'(  Xa)  -  iw 
x  W(Xa)  -  v  W'(  Xa)[l  -f-^'(.v)],  and  g(x)  is  the 
solution  to 

g(x)+v£tH'(Xa)ll  +g<x)1=2aF(Xa),  (10) 

subject  to  g(x  +  2ir)  =g(x).  (The  electric  field  E  is 
given  by  -  ( 2n )  ~ '  f'wdx  g(  x ) .] 

Equations  (9)  and  (10)  were  solved  numerically  and 
the  results  used  with  (5),  (7),  and  <8)  to  determine  hf 
above  threshold.  Below  threshold,  Eq.  (2)  was  in¬ 
tegrated  numerically  to  a  stationary  state  and  then 
linearized  to  calculate  bf  to  O(^)  in  (7).  Since  the 
relative  shift  in  resonant  frequency  is  small,  the  real 
and  imaginary  parts  of  bf  directly  give  and  Q~\ 
respectively,  each  to  within  a  multiplicative  constant. 

These  quantities  are  shown  in  Fig.  1  as  a  function  of 
voltage,  or  £  W  was  chosen  to  be  a  Gaussian  of 
height  1  and  width  6  and  F  had  a  maximum  value  of 
2.2  and  resulted  from  a  repulsive  Gaussian  potential  ot 
width  3,  with  q  =*  777  and  <0  =  0. 1 2. 


0  2  4  6  8  iO 


E/E, 

FIG.  1.  Electromechanical  properties  of  the  incommensu¬ 
rate  lattice-CDW  model:  relative  shift  in  the  Young's 
modulus  ( A<uii/i.ifl)  and  internal  friction  1  Q~'  >.  boih  plotted 
in  arbitrary  linear  units 
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FIG.  2.  Electromechanical  properties  of  TaS?  at  96  K 
(Ref.  12)  The  vertical  bars  indicate  the  maximum  scatter  in 
the  data  points. 


Experimental  results12  are  shown  in  Fig.  2.  The  cal¬ 
culated  Aoj<5  reproduces  the  general  form  of  the  data. 
It  also  contains  a  peak,  however,  which  is  not  apparent 
in  these  data,  although  much  smaller  peaks  have  now 
been  seen  in  some  samples.12, 13  This  difference 
between  theory  and  experiment  may  be  due,  for  exam¬ 
ple,  to  the  existence  of  a  phason  gap  in  the  mean-field 
theory  (see  below),  to  the  one-dimensionality  of  the 
model,  or  to  inhomogeneities  in  the  crystal. 

The  rapid  rise  of  Q~{  above  threshold  is  clearly 
seen,  with  Q~'  reaching  a  maximum  well  before  u»5 
has  leveled  off,  a  feature  emphasized  in  the  original 


experimental  report.2  The  gradual  decline  of  Q  1  at 


higher  voltages  was  not  reported  prior  to  these  calcula¬ 
tions,  and  has  by  now  been  repeatedly  observed.12, 13 

A  qualitative  interpretation  of  the  Young’s-modulus 
results  is  that,  for  small  q,  the  rigidity  of  a  CDW 
pinned  to  a  lattice  enhances  the  rigidity  of  the  lattice. 
For  fields  larger  than  threshold,  this  mechanical  cou¬ 
pling  of  the  two  systems,14, 15  and  hence  the  enhance¬ 
ment  of  the  Young’s  modulus,  diminishes. 

This  decoupling  of  the  sliding  CDW  from  the  lattice 
excitations  can  be  demonstrated  particularly  dramati¬ 
cally  with  the  choice  of  infinite-range  CDW  internal 
interactions.  In  this  case  the  CDW’s  bare  excitation 
(phason)  spectrum  has  a  gap:  wphason  ■/*  0  as  <7  —  0. 
When  the  CDW  is  pinned  to  the  lattice  this  results  in  a 
gap  in  the  spectrum  of  the  combined  lattice-CDW  sys¬ 
tem.  [This  was  checked  numerically  by  determining 
that  nfiui.q  — ■ •  0)  -/•  0  below  threshold.!  Above 
threshold,  comparison  of  (9)  to  the  negative  of  the 
derivative  of  (10)  with  respect  to  .<r  shows  that 
h(x)—  -tf  (.v)  as  q.w—  0.  Use  of  (4),  (5),  (7), 
and  (8)  then  gives 


6  /  (  </,  CM  —  0 ) 


f  *1 

J  -  IT 


dx 


2tt«  -  v  -  g ( .v ) )  =  0. 
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since  g  is  periodic  and  continuous  above  threshold. 
Thus  the  rigidity  of  the  CDW  is  coupled  strongly  to 
the  lattice  below  threshold,  in  this  case  producing  a 
gap,  but  less  strongly  above  threshold,  as  indicated 
here  by  the  disappearance  of  the  gap.16 

When  we  turn  to  the  internal  friction,  below  thresh¬ 
old  many  IDF's  are  suppressed  by  the  pinning,  whiie 
above  threshold  more  can  be  excited.  These  dissipate 
energy  and  enhance  Q~l.  As  the  CDW  velocity  in¬ 
creases  further,  however,  the  coupling  between  the 
lattice  and  the  CDW  weakens  and  the  excitation  of 
IDF's,  and  hence  Q~ *,  decreases  slowly. 

The  importance  of  CDW  IDF’s  can  also  be  demon¬ 
strated  by  our  leaving  them  out,  for  example,  by  mak¬ 
ing  the  CDW  stiffness  coefficients  Df  in  (2)  very 
large,  so  that  (Jj-{U)  and  Uj  vanish.  One  then  re¬ 
covers  a  model  similar  to  that  studied17  to  illustrate  a 
symmetry-breaking  effect  of  the  CDW  current.  From 
(7)  and  (8K  8/  in  this  limit  is,  to  0(*z), 
-  / ct»  ( 2  7r ) —  1 J* '  dx  2^  W(  2wa  -  x ) ,  a  constant.  Thus 
to  Oil u.)  a  rigid  CDW  produces  no  voltage  dependence 
in  the  crystal’s  mechanical  properties;  hence,  the  ef¬ 
fects  calculated  in  Ref.  17  were  0(/u2).  The  presence 
of  IDF’s,  however,  allows  a  first-order  contribution,  as 
shown  in  Fig.  1.  Since  the  crystal  lattice  is  much 
stiffer  than  the  CDW,  first-order  effects  are  much 
larger  than  second-order  effects,  and  IDF's  of  the 
CDW  thus  produce  voltage-dependent  shifts  in 
mechanical  properties  of  magnitude  larger  than  those 
predicted  with  no  IDF's. 

The  dependence  of  lattice  properties  on  the  CDW 
velocity,  v,  is  also  of  interest.  Without  IDF’s  the 
sound  velocity,  s,  is  found17  to  be  an  analytic  function 
of  v  at  v  —  0,  with  an  antisymmetric,  linear  leading 
term  and  higher-order  corrections  that  are  much  small¬ 
er  if,  as  is  the  case  in  CDW  experiments,  u  «  5.  The 
exhancement  of  the  Young’s  modulus  studied  in  this 
Letter,  however,  will  produce  a  symmetric  in¬ 
dependent  component18  in  j(u),  which  should  there¬ 
fore  deviate  from  antisymmetry  over  observable  CDW 
velocity  scales  much  less  than  s.  The  inclusion  of 
IDF's  also  results  in  the  v  —  0  limit  being  a  singular 
limit,  the  threshold  at  which,  for  example,  the  electric 
field  is  related  to  v  in  a  nonanalytic  fashion.5  and  the 
function  g  becomes  singular.  It  is  thus  possible  that  8/ 
in  (7),  and  hence  the  sound  velocity,  will  also  be  non¬ 
analytic  functions  of  v.  Thus,  in  addition  to  strongly 
enhancing  the  current- induced  shifts  in  lattice  proper¬ 
ties,  the  CDW  IDF’s  will  contribute  a  significant  sym¬ 
metric  component  to.  and  possibly  also  change  the 
analyticity  of,  the  sound  velocity  as  a  function  of  CDW 
current. 

A  single-domain,  incommensurate,  classical,  man\- 
body  model  of  the  type  that  has  exhibited  a  variety  nt 
electrical  properties  of  CDW  conductors4-"  has  been 
solved  and  seen  to  account  for  electromechanical  prop- 
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ernes  as  well  As  is  the  ease  for  the  electrical  proper¬ 
ties.  the  internal  degrees  of  freedom1415  of  the  CDW 
were  seen  to  play  an  essential  role. 
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mA  /-lO/Sr)N)  >(  Xj) 

-/uY(2n-A/rl(3/art f  "  ytx  hix  =* 0. 
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We  report  the  discovery  of  an  oscillatory  instability  in  the  dynamics  of  incommensurate  structures. 
The  isciilations  survive  the  thermody  namic  limn.  The  instability  occurs  for  both  lone-  and  short-range 
interactions  The  frequence  and  stability  of  the  oscillations  are  studied. 
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Quasiperiodic  or  incommensurate  systems  have  at¬ 
tracted  considerable  interest  in  a  wide  variety  of  con¬ 
tests,  including  adsorbed  layers,  structures  of  solids,  the 
onset  of  chaos,  and  localization  (see,  for  example.  Refs. 
1-3).  They  also  exhibit4  a  variety  of  the  nonlinear  prop¬ 
erties  of  sliding  charge-density  waves  (CDW’s),  includ¬ 
ing  dc  characteristics,  ac  response,  ac-dc  interference, 
ind  electromechanical  properties. 

We  report  the  discovery  of  a  new  phenomenon,  a  bulk 
oscillatory  instability,  in  the  dynamics  of  sliding  incom¬ 
mensurate  structures:  and  we  study  some  of  its  proper- 
ties  The  existence  of  bulk  oscillations  is  surprising  since 
the  general  belief  has  been  that  the  phase  of  any  oscilla¬ 
tions  in  a  sliding  noncommensurate  system  would  vary 
through  the  sample,  thus  cancelling  the  oscillations  in 
the  thermodynamic  limn,  as  indeed  occurs  in  perturba¬ 
tion  theory.5''  The  instability  is  therefore  breaking  the 
translational  symmetry  of  the  bulk  noncommensurate 
system.  Finally,  in  the  light  of  this  new'  result,  we  dis¬ 
cuss  the  long-standing  problem  of  oscillatory  voltage 
tluctuations  in  CDW  conductors.7'" 

The  first  model  we  consider  is  one  for  whose  dc  prop¬ 
erties  an  exact  solution  is  available,4  namely,  a  set  of  .V 
particles,  all  interacting  equally  with  each  other,  subject 
to  a  pinning  force  P Lv )  and  a  uniform  force  F. 

L'j  -P<Hj  +  L))+‘L  -l',+F.  (1) 

Here  //  is  the  lattice  spacing;  L',  is  the  displacement  of 
the  yth  particle  (so  that  Hj  +  Uj  is  its  position); 
L  ”.V  "  P<x  +  2rr)  wmP(x );  and  we  study  the 

limit  .V  -  oo  an(j  Hih 7  “an  irrational.  Because  CDW's 
are  overdamped.'  "  purely  rclaxationai  dynamics  is 
used.  This  model  can  be  thought  of  as  a  mean-field 
theory  for  the  incommensurate  chain.  Following  a  sim¬ 
ple  construction  ! '  it  can  be  shown,  however,  that  it  can 
1 1  so  be  considered  a  mean-livid  theory  of  a  CDW  subject 
to  a  spatiallv  random  distribution  of  identical  pinning 
venters  (With  in  finite -range  interactions,  the  distinction 
between  the  quasiperiodie  phase  variable  Hj  and  a  com¬ 


pletely  random  phase  variable  pj  is  irrelevant.)  It  has 
also  provided4  accounts  of  ac-dc  interference  experi¬ 
ments  in  the  CDW'  compound  TaSi,  and  the  scaling  of 
field-  and  frequency-dependent  conductivities. 

The  exact  dc  solution  of  this  model  is4  (/)=*[•/ 
+  g(Hj  +  rt)  where  g(.v)  is  the  solution  of  the  bound¬ 
ary-value  problem 

c(l  +g')  =/’(.v  +  g)  —  g,  g(x  +  2rt)  — g(. v ). 

This  solution  was  used  here  as  the  initial  configuration 
for  a  numerical  integration  of  (1).  with  ,V  particles,  and 
it  was  checked  that  the  numerical  procedure  was  stable. 
By  addition  of  any  small  perturbation  the  dynamic  sta¬ 
bility  of  this  solution  can  be  studied.  The  voltage  (F) 
versus  time  plot  in  Fig.  1  shows  the  result  of  such  a 
study.  The  results  presented  in  Figs.  1  and  2  are14  for 
P  (x)  =”8sinx  +  1 2sin4.x. 

It  is  immediately  apparent  that  the  dc  solution  is 
dynamically  unstable,  with  oscillatory  fluctuations  grow¬ 
ing  exponentially  at  first  and  then  saturating.  The  ap- 


F1G.  I  Voltage  vs  lime  for  mean-field  theory,  with  r"“IO. 
Note  the  exponential  divergence  followed  by  saturation, 
characteristic  ol  an  instahilits  in  a  nonlinear  system 
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CDW  CURRENT („) 

FIG.  2.  Frequency  vs  velocity  for  the  Frenkel-Kontorova 
model.  The  segments  satisfy  i7ci>”(///2jt)*  where  n—  2  for 
the  largest  segment  and  cascades  through  n  ~  1 ,2.3.4  as  the  ve- 
locity  is  decreased  further 

pearance  of  harmonic  content  in  Fit  ),  as  well  as  satura¬ 
tion.  are  both  due  to  nonlinearity  and  they  are  seen  to 
occur  at  the  same  time,  as  expected.  Time  series  ob¬ 
tained  with  ,V  =■  1 44,  233,  and  377  are  essentially  indis¬ 
tinguishable.  Thus  the  results  shown  in  Fig.  1  describe 
the  thermodynamic  limit  and  are  not  a  finite-size  effect. 
Moreover,  a  finite-size  effect  would  not  be  expected  to 
show  the  exponential  divergence  seen  in  Fig.  1,  which  is 
instead  the  signature  of  a  dynamic  instability. 

In  addition,  we  performed  a  stability  analysis  of  the  dc 
solution  to  (1)  in  the  limit  .V— *  °°.  CDW  experiments 
are  often  current  driven,  and  so  we  considered  (1)  in  the 
presence  of  "normal  electrons"  by  keeping  fixed  a  total 
current,  i-  +  onF.  where  an  is  the  conductivity  of  the 
linear,  normal  channel.  It  can  be  shown  that,  for  any 
pinning  potential,  oscillations  at  (complex)  frequency  o> 
will  occur  under  conditions  of  fixed  total  current  v  +  o„F 
only  if  a>  is  a  root  of  —  tyn  where  cr(a>)  is  the 

linear  response  function  of  the  sliding  structure.  Thus, 
while  an  instability  under  conditions  of  fixed  voltage  F 
would  correspond  to  a  pole  of  a  crossing  into  the  upper- 
halt  complex  o)  plane,  the  instability  shown  in  Fig.  1.  at 
fixed  current  r,  corresponds  to  a  zero  of  cr  crossing  into 
the  upper-half  plane.  We  solved  the  linear-response 
equation  and  obtained  the  zeros  of  <r  for  different  veloci¬ 
ties.  At  r="!0  the  result  w-as  a  normal-mode  frequency 
whose  real  and  imaginary  parts  both  agreed  precisely 
with  the  diverging  oscillation  seen  in  Fig.  I,  confirming 
that  it  shows  an  oscillatory  dynamic  instability  inherent 
in  the  mean-field  theory  of  the  dynamics  of  incommensu¬ 
rate  structures. 

V>e  calculated  the  complex  normal-mode  frequency  at 
different  velocities.  As  the  velocity  decreases  below  II, 
the  zero  moves  into  the  upper-half  complex  plane,  signal- 
tng  in  instability  of  bifurcation.  There  is  thus  a  critical 
velocity  r,  above  which  the  dc  solution  is  stable  (oscilla¬ 
tions  decay)  and  below  which  the  oscillations  persist. 

I  oi)4 


The  instability  was  observed  for  a  wide  range  of  on,  with 
t\  decreasing  as  <r„  increases. 15 

The  solution  of  the  linear  fluctuation  equation  at  r,., 
Si ',i  =  q(Hj  +  vt  )e  provides  some  rudimentary  in¬ 
sight  into  the  origin  of  the  instability.  The  exact  solution 
of  the  dc  motion4  showed  that  there  are  values  of  the 
pinning  potential  (essentially  the  peak  values)  which  a 
locally  stable  static  solution  avoids  by  having  discon¬ 
tinuities  in  the  function  g < a  )  defined  above.  Static  solu¬ 
tions  do  exist  with  particles  in  these  regions,  but  such 
solutions  are  unstable.4  A  sliding  system,  however,  must 
have  particles  in  these  regions  and  a  continuous  g(x), 
with  dg/dx  sharply  peaked  in  these  regions  for  small  e. 
We  find  that  the  unstable  fluctuation  q(x)  has  |  rj  |  larg¬ 
est  just  where  dg/dx  is  peaked.  Thus  the  oscillatory  in¬ 
stability  of  the  dc  solution  may  be  in  some  measure  a  dy¬ 
namic  consequence  of  the  existence  of  unstable  static 
states. 

To  see  whether  the  oscillatory  instability  exists  in 
finitely  coordinated  systems  we  studied  a  minimally  co¬ 
ordinated  system:  an  incommensurate  chain  with  only 
nearest-neighbor  interactions,1 

i'j  2i'j  +  + 1  +  F.  (2) 

The  incommensurate  chain  has  been  shown4  to  give  an 
account  of  the  difference  between  the  ac  and  dc  interfer¬ 
ence  properties  of  NbSej  and  TaS3.  A  previously  ob¬ 
tained4  dc  solution  to  (2)  was  .used  as  an  initial 
configuration  in  a  numerical  integration  with  .V  particles. 
Increasing  .V  sufficiently  produced  identical  plots,  ensur¬ 
ing  that  the  results  represented  the  thermodynamic  limit. 
Different  pinning  potentials  were  studied.  So  long  as 
P(x)  produced  a  threshold,  whether  or  not  P(x)  con¬ 
tained  harmonics,  resulting  voltage-time  plots  showed 
precisely  the  exponential  instability,  followed  by  satura¬ 
tion.  that  was  seen  in  Fig.  I  for  the  mean-field  case. 

Thus  the  instability  occurs  at  coordination  number 
infinity  and  two.  It  is  therefore  expected  to  occur  in  in- 
commensuratelv  pinned  systems  at  all  intermediate  coor¬ 
dinations  in  one.  two,  and  three  dimensions. 

The  oscillation  frequency  as  a  function  of  v.  for 
nearest-neighbor  interactions  and  the  same  P(x)  as  in 
Fig.  1.  is  shown  in  Fig.  2  and  has  some  interesting 
features.  It  is  predominantly  linear.  Further,  although 
the  large-r  limit  has  transients  with  the  trivial  frequency 
<u“r  which  would  be  exhibited  by  a  single  particle  in 
P f.v).  once  the  instability  takes  over  (c<9)  a  new 
characteristic  frequency  and  length  appear  which  are 
determined  not  only  by  Pix),  but  also  by  the  sliding 
structure  itself.  In  the  large  linear  region  in  Fig.  2. 
<7<o=*(/f/2.T)-.  That  is.  the  oscillatory  instability  re¬ 
veals  the  length  scale  H  provided  by  the  lattice  spacing, 
which  is  the  length  scale  corresponding  to  the  wave¬ 
length  of  a  CDW 

Further,  as  r  is  decreased,  there  are  a  number  of  first- 
order  transitions,  where  r/w  changes  abruptly,  being 
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glven  hs  <//-2.r)".  n=l.  1.2.3.  .  as  the  velocity  is  de¬ 
creased. 

The  properties  of  the  Erenkel-Kontorova  model  [Hq. 
12)1  have  attracted  a  great  deal  of  attention  (see.  e.g.. 
Rci  2).  L  is  no«  clear  that  its  dynamic  properties  are 
considerably  richer  than  previously  realized. 

Turning  to  CD  A  conductors,  the  dominant  source  of 
pinmniz  is  believed  to  be  a  random  potential  due  to  lat¬ 
tice  detects.  While  Eq.  •  !)  is  also  a  mean-field  theory  of 
random  pinning,  we  do  not  yet  know  whether  a  bulk  os¬ 
cillators  instability  also  occurs  in  systems  with  random 
pinmnsi  and  short-range  interactions. 

No  critical  velocity .  iv.  has  been  reported  in  the  exper¬ 
imental  literature.  Some  experiments  in  the  time 
domain  have,  however,  revealed  transient  oscillations.  6 
While  the  existence  of  transients  has  no  clear  explana¬ 
tion  within  the  finite-size 1  '■ 1  or  contact  "s  theories,  a  nat¬ 
ural  interpretation  in  the  present  context  is  that  in  these 
experiments  r  >  iy  so  that  fluctuations  oscillate  but  de¬ 
cay  . 

The  origin  of  CDW  voltage  oscillations  remains  an 
open  question.  The  unexpected  observation  of  bulk  oscil¬ 
lations  in  models  which  have  exhibited  a  wide  variety  of 
other  properties  of  sliding  CDW’s  may  be  the  basis  for 
the  answer,  or  merely  a  tantalizing  coincidence. 

In  conclusion,  we  have  reported  the  discovery  of  new 
phenomenon — an  oscillatory  instability — in  the  dynam¬ 
ics  of  incommensurate  structures. 
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Mixinn  in  chars>e-density-wave  conductors 
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Harmonic  and  direct  ac  mixing  properties  <>l  the  Fukuvama-I.ee-Ricc  model  and  the  incom¬ 
mensurate  chain  are  determined  and  compared  with  experimental  data.  Field  and  frequency 
dependences  ol  both  the  amplitudes  and  phases  ot  both  ol  these  responses  are  examined  Agree¬ 
ment  with  experiment  is  generally  good.  For  example,  classical  models  can  vieid  low  harmonic- 
mmng  quadrature  components  simultaneously  with  substantial  frequency  dependence  ;n  both 
components  of  the  linear  response. 


An  interesting  series  of  experiments,  involving  the  non¬ 
linear  mixing  of  ac  signals  in  charge-density-wave 
<CDW)  conductors,  was  initiated  by  Seeger.  Maver.  and 
Philipp •  and  thoroughly  extended  by  Miller  and  co- 
workers.'"'  These  experiments  merit  theoretical  study  in 
their  own  right  as  probes  of  the  unique  properties  of  CDW 
conductors. 

F:urther,  the  suggestion  has  been  made  that  the  results 
"may  prove  difficult  to  reconcile  with  any  classical 
theory."4  It  is  thus  important  to  determine  the  ac  mixing 
properties  of  classical  models  of  sliding  CDW's  to  see 
w  hether  there  is  indeed  a  failure  of  the  classical  picture  of 
bulk  CDW  motion,  and  whether  these  experiments  give  us 
evidence  for  Bardeen's  fascinating  proposal6  that  CDW 
conductors  are  exhibiting  macroscopic  quantum  tunnel- 
ing. 

The  principal  experiments2*5  are  direct  mixing  <or 
rectffication)  and  harmonic  mixing.  The  sample  is  driven 
by  a  voltage  of  the  form 

E  (r )  ”£o  +  Ei  cos(a>i/  +<?)  +  E  ^cosfwsr )  ,  0) 

where  E\  and  £N  are  generally  small  compared  to  the 
threshold  voltage  £7-.  The  two  measurements  consist  of 
detecting  the  component  of  the  current  with  frequency  ton. 
where  \  ~  co;  for  direct  mixing,  and  oin  =  2cui  —  w; 

for  harmonic  mixing.  The  frequencies  are  chosen  so  that 
cut)  is  much  smaller  than  ai|  and  tur 

The  properties  of  two  different  models  are  reported 
here:  the  Fukuyama-Lee-Rice  fFLR)  model  of  random 
pinning  and  the  incommensurate  chain.5  Perturbation  1 
theory  was  used  to  obtain  analytic  information  for  the 
FI.R  model  at  large  bias  fields  E q.  The  incommensurate 
chain  was  simulated  numerically  to  obtain  solutions  in  the 
strong-coupling  region  closer  to  the  sliding  threshold  field 
£/.  The  equation  of  motion  for  the  incommensurate 
chain  with  infinite-range  interactions1011  is 

—  —iv  -u,  +  Psm(Hj  +  u,  )  +  £(/)  .  <2) 

Lit 

where  it,  is  the  displacement  of  the  yth  particle,  u.  is  the 
center-of-mass  displacement.  P  is  the  pinning  strength, 
and  I/!<2.r)  is  cnosen  to  be  f  V 5  —  1  1/2. 

By  vimuiatinu  Fq.  <21  in  systems  up  to  277  particles  in 
size  we  were  able  to  restrict  linite-M/e  effiects  to  the  region 
vers  close  to  threshold  and  ensure  that  'he  results  ob¬ 
tained  at  dl  other  fieids  reliably  retlected  'he  thermo- 

<s 


dynamic  limit.  The  applied  frequencies  were  chosen  to 
have  a  rational  ratio  so  that  the  response  (in  the  thermo¬ 
dynamic  limit!  is  periodic.  Transients  were  allowed  to  de¬ 
cay  for  a  time  to =  T  sufficiently  long  that  to^T  and 
tt)  =  2T  give  results  that  are  indistinguishable  within  nu¬ 
merical  error.  The  sampling  rate  of  the  time  series  yvas 
chosen  high  enough  that  increasing  it  further  produced  no 
noticeable  changes. 

We  now  compare  the  theoretical  results  with  those  of 
experiments,  examining  each  feature  of  the  data  in  turn. 

1.  HARMONIC  MIXING  PHASE  SHIFT 

The  experimental  result  which  has  been  most  em¬ 
phasized"'4  ■  is  the  failure  of  some  experiments  to  observe 
an  "internal"  phase  shift  in  the  small  difference-frequency 
harmonic-mixing  response.  Classical  models  were  conjec¬ 
tured"  to  yield  nonzero  phase  shifts  at  high  applied  fre¬ 
quencies  tot.  Wonneberger12  then  showed  that  the  classi¬ 
cal  single-particle  model  exhibits  a  zero  phase  shift  for  all 
<o i  at  large  dc  bias  fields.  Fie  proved  this  result  to  leading 
order  in  perturbation  theory  and  argued  in  an  appendix 
that  the  result  is  valid  everywhere  above  threshold. 
Nonetheless,  it  was  claimed  that  the  nonobservation  of  a 
phase  shift  above  threshold,  where  in  addition  both  com¬ 
ponents  of  the  linear  ac  response  have  substantial  frequen¬ 
cy  dependence,  may  be  difficult  to  reconcile  wuh  any  clas¬ 
sical  theory  4  It  was  further  claimed  that  this  experimen¬ 
tal  observation  provided  particularly  significant  evidence 
for  CDW  tunneiimt  " 

We  determined  the  harmonic-mixing  properties  of  the 
FLR  model  using  perturbation  theory  at  large  dc  bias 
fields.  The  principal  resuit  is"'  that  the  response  com¬ 
ponent  at  trequency  o™  is  proportional  to  cosIom  4-2 ol  in 
the  limit  of  small  con  for  large  dc  bias  fields  /.' , >  and  all  on. 
The  sine  component  vanishes  iincariy  with  um.  Thus  the 
classical  FLR  model  also  exhibits  t  zero  internal  phase 
shift  at  large  bias  fields. 

In  that  region,  however,  the  linear  ac  response  has  no 
strong  Irequency  dependence.  It  is  therefore  necessary  to 
prone  the  region  of  lower  vie  bias  closer  to  threshold.  For 
this  purpose,  we  turn  to  the  results  tor  the  incommensu¬ 
rate  chain 

Figure  I  shows  some  typical  held  and  .a,,  dependences 
oi  both  components  ol  the  harmonic-inixine  -ien.il.  It  ,s 
seen  that  the  un-oi-phase  component  :s  considerably 
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for  ’be  incommensurate  chain.  Figure  2(h)  snows  the  ex¬ 
perimental  data.-'  ana  Fig.  2(c)  -.hows  the  resuits-4 
a  men  ha\e  been  presented  based  on  quantum  tunneling 

arguments. 

'..v  The  peak  \  In  Fig.  2(a).  the  positions  (fields)  and 
:r.e  reients  ot  the  peaks  change  in  the  same  sense  with  fre¬ 
quency  as  seen  experimentally.  The  present  resuits  lor  the 
incommensurate  chain  change  somewhat  more  rapidly 
with  frequency  than  do  the  experimental  data. 

'hJ  The  threshold.  The  incommensurate  chain  is  seen 
in  Fics.  2(a)  and  2(b)  to  give  the  experimentally  observed 
qualitative  forms  in  the  vicinity  of  threshold.  The 
threshold-lowering  effect  of  an  appiied  ac  field  is  also  ap¬ 
parent  in  the  incommensurate  chain  results,  the  effect  in¬ 
creasing  as  the  ac  frequency  decreases,  as  observed  experi¬ 
mentally  .  ^ 

A  comparison  with  Fig.  2(c)  then  shows  that  these 
threshold  features  reveal  substantial  differences  between 
the  properties  of  the  classical  chain  and  the  results  of  the 
quantum  tunneling  analysis. 

<c)  The  negative  dips.  The  experimental  data  [see  Fig. 
2 f b ) ]  mow  that  for  lower  frequencies  u>i  the  harmonic¬ 
mixing  response  goes  negative  above  a  (frequency- 
dependent)  field  (the  phase  switches  rapidly  from  a  value 
at  or  near  zero  to  one  at  or  near  180°).  Thereafter  the 
response  approaches  zero  from  below.  Above  a  certain 
frequency,  however,  this  no  longer  happens  and  the 
harmonic-mixing  response  is  always  positive. 

This  points  up  another  qualitative  difference  between 
the  two  sets  of  theoretical  results:  While  the  incommensu¬ 
rate  chain,  like  the  experimental  results,  shows  an  « t 
move  which  the  response  no  longer  changes  sign,  the 
quantum  tunneling  analysis  appears  to  show  a  sign  change 
for  all  wi. 

f l)  Below  threshold.  The  calculated  amplitude  of  the 
harmonic-mixing  signal  is  about  two  or  three  orders  of 
magnitude  less  than  the  peak  above  threshold,  which  is 
consistent  with  the  experimental  data  [see  Fig.  2(b)], 


0.025 -  7  -EC'  ‘  -3 


0020  -i 


0.000  -I - I - r- - . - ■ - - - - - — - . 

0.0  05  V0  ’5  2.0  25  3.0  35 

oc  vein  Ec 


RECTIFICATION 


vbi„!mv> 


FIG.  3.  dc  field  (£o)  and  frequency  Uu i )  dependencies  of  the 
direct-mixing  current  (,V  ' 1  V  Ju,/dt )  magnitude:  <a)  the  in¬ 
commensurate  chain.  .V— 233.  P  —  3.0.  and  I:  ■  =*£;  =0.30;  (h) 
experiment  ( Ref.  2). 


penmental  data  are  seen  to  be  reasonably  well  reproduced 
bv  the  incommensurate  chain. 


III.  DIRECT-MIXING  PHASE 

Experimentally,-"*  the  response  at  —  <hp  to  the 

cosine  input  signal  of  Eq.  ( 1 )  is  found  to  he  proportional 
to  cosGoii/  +o)  lor  small  u> <>.  A  leading-order  perturbative 
dilution  of  the  Fi  R  model11  showed  this  same  feature. 
This  conclusion  can  be  seen  to  be  true  to  all  orders  as  fol¬ 
lows  Consider  the  case  o=0.  Equation  ( 1 )  is  then  sym¬ 
metric  with  respect  to  interchanging  iO)  and  ox  The  total 
•esponse  current  at  frequency  plus  or  minus  aiu.Jnsh  then 
iiv>  :ias  this  symmetry.  Writing  /dm  W,  siniwn/ ) 
-/  ci is Iw../ 1  'hen  requires  *  —•/,(«;,  u>i ). 

■fat  I  *0  as  <o i  -  a-  Even  if  the  input  tields  ire 
r.muea  to  sines  n  shiftmu  the  origin  of  time,  the  !ow- 
‘reaven.’.  direct -mixirtL’  response  is  stiil  cosine. 

I\.  OIRl  (  T- MIXING  M  M.NI  H  1)1 

i  :  ’  ,re  :,.i)  .ro,  and  frequent',  dependence  I 
'!■■  a  iL'iutuile  ■!  nreet-iimim.'  espouse  tor  'he  ,u- 

■  oe  or.  I  i.rure  >(bi  nows  ■  (-.e  , or respond - 

■'  i'eriuiei.'ai  t  Put  orinar  o  lie 


V.  C  OM  El  SION 

The  classical  incommensurate  chain  has  been  seen  to 
provide  a  fairly  complete  account  of  the  heid  and  treuuen- 
cy  dependencies  of  both  amplitude  and  phase  components 
of  the  direct  and  harmomc-mixinu  responses  of  sliding 
CDW's. 

For  example,  the  challenge  to  classical  theories  posed" 
by  some  nonobservations  of  the  ••.armomc-mixine  phase 
shift  has  been  resolved  by  the  observations  that  1  a  l  phase 
shuts  nave  been  seen  experimental!'.,  ami  'hi  classical 
theories  can  account  for  phase  shift'  as  'tnui!  as  the  exper¬ 
imental  uncertainties.  ’It  :n,iv  be  interesting  now  \o  have 
both  exrerimeni.ii  ana  theoretical  pnase  smtts  determined 
more  prcciseis  to  provide  more  teiiing  comparisons  >1 
T.eories  and  experiment 

I  he  pinnine  •(  <  I)A  s  •x-iicvi:  ■  ■  'w  tue  '  ■  :  imiom 

mp-inties.  me.  i. >r  i  periodic  nr.ciiMir.i','  poteii- 

■ ,  !"  end  ui  i  question.  ■  pen.  m-  c  'c  '".at  '::c  issue  a 

vie.  Tiv.'inme.’tsurate  vsiem.s  repr- since  .kp  :  inetv  «t 
a'  erv -'il  <  |)W  •-p.ivr’ie.  s  cm  '.''  c  ■  per  v  macro.  ,Hj 
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smaller  than  the  ;n-phase  component,  but  appears  to  he 
approaching  a  nonzero  limit  as  con  approaches  zero.  \s  om 
decreases  through  0.05.  0.025.  and  0.0 1 25  the  >ut-oi- 
phase  component  shows  practically  no  change  on  the  scale 
of  the  plot. 

This  limit  is  generallv  small.  Kor  the  data  shown  in  Fig. 
1.  (oiya>\  is  about  0.01.  it  corresponds  to  a  phase  shift  of 
about  10°.  which  is  the  reported  resolution  of  the  experi¬ 
ments.45  Experimentally  luqj'id]  was  chosen  to  be  about 
0.001.  We  found  that  this  phase  shift  decreases  as  £i  and 
E :  are  increased  with  the  other  parameter  fixed.  The 
value  of  oji  in  Fig.  1  was  chosen  approximately  to  maxim¬ 
ize  the  phase  shift. 

Earlier  work11  has  established  substantial  frequency 
dependence  in  both  components  of  the  linear  ac  response 
in  this  region.  One  may  thus  conclude  that  classical 
theories  can  indeed  account  for  the  absence  of  strong 
quadrature  in  the  harmonic  mixing  and  simultaneous 
presence  of  strong  frequency  dependence  in  the  linear 
response. 14 

The  above  classical  results  do  appear  to  show  a  nonzero 
quadrature,  however,  and  if  a  classical  picture  were  ap¬ 
propriate,  one  would  expect  to  see  a  phase  shift  under 
some  experimental  conditions.  A  harmonic-mixing  phase 
shift  has,  in  fact,  been  reported. 15  The  ratios  E\/Et  and 
EilET  were  smaller  than  in  the  experiments2,3  that  saw 
no  phase  shift.  This  is  consistent  with  our  observation 
above  that  the  phase  shift  decreases  as  E\  and  Ei  in¬ 
crease. 

Turning  to  the  region  below  threshold,  the  calculated 
out-of-phase  component  is  found  to  be  much  larger  than 
the  in-phase  one  for  small  tuo-  This  is  consistent  with  ex¬ 
perimental  reports4  of  an  internal  phase  shift  of  about 
n/2. 


HARMONIC  MIXING 
INCOMMENSURATE  CHAIN 


m-phase 


6c  field  E0 

EIG  I  uc  tieid  1 E; :  ind  difference-frequency  Uan)  depen¬ 
dencies  of  both  components  of  : he  harmonic-mixing  response 
current  f.V " 1  ol  the  incommensurate  chain 

=*10.  /’  “  *0.  and  /.';  -  /:  .  “0.2 )  \  1 1  quantities  are  dimension¬ 

less  [see  Eq  <2)1.  The  curves  lor  V  “  »7~  indicate  the  a/e  a 
truncation  error  Both  components  have  ’mite  :>n  -O  limits, 
with  i he  out-phase  espouse  considerable  .mailer  than  he  n- 
ph.ise  response 


IE  H  \R MON  1C -MIXING  MAGNITl  I)L 

Since  the  harmonic-mixing  phase  shift  is  generally 
small,  we  locus  now  on  the  magnitude.  Figure  2(a)  shows 
the  field  'En)  and  frequency  (w\)  dependencies  of  the 
magnitude  of  the  harmonic-mixing  response,  as  calculated 
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ABSTRACT 

A  substantial  disagreement  between  mode-locking  experiments,  and 
predictions  based  on  a  classical  model,  were  previously 
interpreted  as  proof  of  the  fundamental  deficiency  of  the 
classical  model.  We  report  classical  calculations  which  disprove 
the  predictions,  and  agree  with  the  experimental  data.  Contrary 
to  the  two  previous  reports,  classical  models  exhibit  complete 
mode  locking  at  high  fields  and  frequencies.  The  importance  of 
nonlinear  resonance  is  stressed.  We  present  results  for  the 
differential  resistance,  spectral  response  and  phase  boundaries. 
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Collective  Nonlinearity  in  Solids 


When  nonlinearity  drives  many  degrees  of  freedom  m  dynamic 
incommensurate  systems,  these  models  reproduce  the  properties  of 
charge  density  wave  conductors  (CDW's)  in  an  unexpectedly  wide 
range  of  experimental  situations.  The  collective  and  nonlinear 
properties  of  these  dynamical  systems  are  thus  frequently  much 
more  important  than  the  precise  nature  of  the  pinning  potential. 
In  this  section,  the  role  of  the  internal  degrees  of  freedom  is 
elucidated  in  a  variety  of  experimental  observations,  for  example 
the  existence  of  the  threshold;  the  shape  of  the  I-V 
characteristic;  ac/dc  coupling  and  interference;  dissipative 
effects;  and  elastic  phenomena. 


a)  The  Threshold. 


The  fundamental  issue  of  the  existence  of  a  threshold  in 
non-commensurate  systems  can  be  thought  of  rather  naturally  in 
terms  of  collectivity.  In  an  incommensurate  system,  a  weak 
pinning  potential  excites  distortions  with  the  wavevector (s)  of 
the  potential,  and  also  their  harmonics.  The  harmonics  decay 
exponentially  in  strength1,  however,  so  that  an  infinite  number 


1  S.  J.  Shenker  and  L.  P.  Kadanoff,  J.  Stat.  Phys.  27, 
631  (1982) 


of  degrees  of  freedom  are  not  excited.  There  is  no  threshold2. 

As  the  pinning  strength  is  increased,  the  regime  of 
exponential  decay  gets  pushed  out  to  higher  and  higher  harmonic 
number  (through  the  nonlinearity) .  At  a  finite  critical  pinning 
strength,  the  harmonics  never  decay  exponentially,  but  only  as  a 
power  law.  At  this  point  all  harmonics  are  contributing  and  the 
problem  is  a  collective  one.  It  is  precisely  at  this  point  that 
the  threshold  becomes  finitel. 

A  random  pinning  potential  excites  all  fourier  components, 
no  matter  how  weak  it  is,  because  it  contains  all  fourier 
components.  Thus  random  pinning  is  a  collective  phenomenon,  no 
matter  how  weak  the  potential.  Correspondingly,  a  random  pinning 
potential  has  a  finite  threshold,  no  matter  how  weak  the 
potential . 

b)  The  I-V  curve. 

We  can  also  see  collectivity  "turn  on"  in  the  I-V  curve. 
Consider  a  straight  line,  S,  through  the  origin  with  slope  equal 
to  the  high  field  differential  conductivity  (see  Fig  5)  .  The 
actual  I-V  curve  is  offset  below  S.  This  "dc  offset"  can  be 
understood  as  due  to  the  excitation  of  a  continuum  of  internal 

2  Peyrard  and  S.  Aubry,  J.  Phys  C  16,  1593  (1983)  and 
references  therein. 


8 

£ 


! 

A 


i' 


degrees  of  freedom  (IDF's).  The  idea  in  the  above  discussion  of 
t  j  threshold  works  here  too.  At  high  fields  the  potential  is 
effectively  weak  (since  the  CDW  is  moving  over  it  too  fast  to 
deform  appreciably) .  The  random  potential  is  always  collective 
and  produces  an  offset  even  at  high  currents. 

The  incommensurate  potential  does  not.  It  can  be  shown  from 
perturbation  theory  that  the  I-V  curve  approaches  the  line  S  at 
high  fields.  That  is,  the  offset  is  not  present  in  low-order 
perturbation  theory.  When  IDF's  are  treated  to  all  orders, 
however'  the  offset  is  seen  to  develop  as  the  field  is  lowered. 
This  is  of  course  also  when  the  nonlinearity  is  exciting  many 
IDF's  and  making  the  problem  truly  collective. 

c)  Interference  Effects. 


The  effect  of  long  range  interactions  on  interference 
effects  provides  two  slightly  more  interesting  examples  of  the 
role  of  collectivity. 

i)  "Universality” 

First,  there  is  an  interesting  comparison  between  the 
properties  of  the  randomly  and  incommensurately  pinned  systems. 
In  the  weakly  nonlinear  regime  (weak  pinning  or  high  fields)  it 
can  be  seen  from  perturbation  theory  that  long  range  interactions 
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How  can  NbSe3  get  around  that?  By  having  strong  short-range 
elastic  interactions.  Short-range  elastic  interactions  favor 
distortions  at  long  wavelengths.  If  the  dominant  distortions  are 
at  wavelengths  much  larger  than  the  CDW  wavelength,  the  washboard 
frequency  will  remain  relatively  well-defined,  and  the  system  can 
behave  non-linearly  and  yet  avoid  the  "catch  22". 

Having  longer-range  interactions,  for  example  Coulomb 
interactions,  however  allows  shorter  wavelength  distortions  to 
play  a  larger  role.  An  infinite-range  interaction  which  is 
independent  of  distance,  for  example,  offers  a  restoring  force 
which  is  independent  of  the  wavelength  of  the  distortion,  and 
thus  does  not  favor  long  wavelength  distortions  at  all.  The 
short  wavelength  distortions  will  then  smear  out  the  microscopic 
frequency,  and  the  interference  features  will  fall  victim  to  the 
"non-linear  catch  22". 

iii)  Experimental  Results 

In  a  similar  fashion  to  the  theory  of  vortex  flow, 
perturbation  theory  was  used  to  predict  interference  features  in 
the  frequency-dependent  conductivities  of  sliding  CDW's3. 
Experimental  confirmation  came  immediately  in  studies  of  NbSe34. 

3  L.  Sneddon  Phys .  Rev.  B29,  725  (1984). 

4  A.  Zettl  and  G.  Gruner  Phys.  Rev.  B29 ,  755  (1984). 


wash  out  the  ac/dc  interference  features  of  randomly  pinned 
systems  (see  below) .  In  the  same  regime,  however,  perturbation 
theory  shows  that  an  incommensurate  chain  will  produce  ac/dc 
interference  features  whether  the  elastic  interactions  are  short- 
or  long-range. 

What  happens  to  incommensurate  chains  in  the  strongly  non¬ 
linear  sliding  regime,  where  the  interference  features  are 
strong?  The  features  become  vulnerable  to  long  range 
interactions  (see  Publication  1)  .  Thus  again,  as  we  approach 
threshold  in  FK  systems  we  enter  a  collective  regime  in  which 
their  properties  imitate  those  of  randomly  pinned  systems. 

ii)  "Non-linear  Catch  22" 

Why  do  long  range  interactions  wash  out  interference 
features  in  randomly  pinned  systems?  The  answer  is  that  the 
interference  is  subject  to  a  "nonlinear  catch  22".  Interference 
is  of  course  a  nonlinear  effect,  and  the  source  of  nonlinearity 
is  the  pinning,  so  for  interference  to  exist  the  system  must  be 
being  distorted  by  the  pinning  potential.  This  distortion, 
however,  locally  and  incoherently  modulates  the  CDW  wavelength, 
and  hence  smears  the  washboard  frequency.  Thus,  to  get  ac/dc 
coupling  at  all,  the  system  must  suffer  a  smearing  out  of  the 
frequency  that  defines  the  interference  feature. 
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The  experimental  results  showed  what  were  called  inductive 
dips:  the  dielectric  constant  can  become  strongly  negative  at  the 
interference  feature.  The  non-perturbative  solution  of  the 
incommensurate  chain  then  showed  us  (Publication  1)  that  these 
inertial  dips, as  well,  could  be  understood  as  the  nonlinear, 
collective  response  of  a  classical,  massless  system. 

TaS3,  however,  did  not  show  these  features4,  and  the 
suggestion  was  made  that  long-  range  Coulomb  interactions  in  the 
semiconducting  TaS3  were  responsible  (Publication  1) .  The 
mechanism  may  be  precisely  the  non-linear  catch  22  outlined 
above.  This  proposal  can  be  probed  a  little  further,  however. 
Coulomb  potentials  will  be  easier  for  the  conduction  electrons  to 
screen  away  if  the  underlying  charge  fluctuations  occur  at  low 
frequencies;  and  also  if  the  normal  conductivity  is  high. 
Mihaly  et  al5  examined  both  the  frequency-  and  temperature- 
dependence  of  the  interference  properties  of  TaS3.  They  found 
that  sharp  features  did  indeed  re-emerge  at  lower  applied 
frequencies  and  also  at  higher  temperatures,  where  the  normal 
conductivity  is  enhanced  in  the  semiconductor.  While  careful 
qualitative  tests  still  remain  to  be  performed,  these  results  are 
clearly  qualitatively  in  agreement  with  the  "catch  22". 
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d)  Viscosity  Enhancement 

Collectivity  also  has  dissipative  consequences.  Because  a 
CDW  carries  condensed  charge,  whenever  it  is  deformed  there  will 
be,  as  noted  above,  charge  accumulation,  which  the  normal 
carriers  will  then  strive  to  screen.  Normal  currents  are 
dissipative,  however,  and  these  normal  screening  currents  were 
predicted6  to  enhance  the  viscosity  or  damping  of  the  CDW  motion, 
the  effect  predicted  to  grow  larger  as  the  temperature  is 

decreased,  roughly  proportionally  to  the  normal  resistivity  in 
the  semiconducting  materials.  The  predicted  effect  on  the 

temperature-dependence  of  the  I-V  characteristics  was  first 
qualitatively  verified  by  Monceau.  About  two  years  after  the 
theoretical  prediction,  another  series  of  experiments  by  Fleming 
et  al.  also  confirmed  that  the  viscosity  had  essentially  the  same 
temperature-dependence  as  the  normal  resistivity. 

e)  Elastic  Phenomena 

The  collective  nature  of  CDW  transport  also  plays  an 

important  role  in  the  elastic  properties  of  CDW  conductors.  A 
rigid  CDW  can  couple  to  the  lattice  modes7,  enabling  a  sliding 
CDW  to  make  the  sound  velocity  anisotropic.  The  effect  would  be 

6  L.  Sneddon  Phys.  Rev.  B29 .  725. 

‘  S.  N.  Coppersmith  and  C.M.  Varma  Phys.  Rev.  B30 ,  3566 
(1984)  . 
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very  small,  and  an  analytic  function  of  the  CDW  current.  When 
the  internal  degrees  of  freedom  of  the  CDW  are  taken  into 
account,  however,  the  effect  becomes  much  larger  (see  Publication 
3).  Mathematically,  the  effect  becomes  first  order  in  the 
lattice-CDW  coupling,  while  without  the  IDF's,  only  second  order 
effects  are  possible.  Physically,  the  IDF's  are  enabling  the  CDW 
to  get  a  "grip"  on  the  lattice,  so  that  to  deform  the  lattice, 
one  must  also  deform  the  CDW,  and  this  contributes  to  the  lattice 
stiffness.  As  the  CDW  starts  to  move  relative  to  the  lattice,  it 
has  less  time  to  for  the  IDF's  to  adjust,  and  the  grip  of  the  CDW 
on  the  lattice  is  not  as  strong,  and  the  stiffness  of  the  lattice 
diminishes,  as  seen  experimentally.  Further,  the  IDF's  produce  a 
sliding  threshold  at  which  the  dynamics  is  singular.  This  is 
consistent  with  the  experimental  observation8  of  a  fractional 
exponent  in  the  Young's  Modulus. 

The  collective  nature  of  the  CDW  transport  also  plays  an 
important  role  in  the  internal  friction  of  elastic  vibrations. 
Below  threshold,  many  of  the  IDF's  are  pinned  by  the  pinning 
potential  and  their  response  to  an  external  driving  is  limited. 
It  is  likely  that  there  is  a  gap  in  their  excitation  spectrum. 
As  the  applied  electric  field  passes  above  threshold,  the  gap 
would  vanish,  and  the  IDF's  suddenly  become  more  easily  excited. 
The  excitation  of  IDF's  dissipates  energy,  and  thus  leads  to  a 
marked  increase  in  the  internal  friction. 
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J.W.  Brill  and  W.  Roark,  Phys  Rev.  Lett.  53,  846  (1984). 


f)  Conclusion 

If  the  function  of  randomness  is  only  to  introduce 
collectivity  and  nonlinearity,  then  we  have  a  rather  satisfactory 
situation.  At  high  fields,  incommensurate  systems  are  not 

collective,  and  so  randomness  is  tssential.  At  high  fields, 
however,  we  can  use  perturbation  theory  to  solve  random  pinning. 
At  low  fields,  incommensurate  systems  become  collective  as  well 
as  nonlinear  and  thus  much  more  realistic.  Further,  it  has 
proved  possible  to  solve  for  many  of  their  properties  in  this 
regime.  Thus  incommensurate  systems  may  well  be,  when  combined 
with  the  HIGH-field  solutions  of  randomly  pinned  systems, 

completing  a  solution  of  sliding  CDW's  by  providing  a  theory 

which  includes  the  most  important  physical  features,  and  yet 
which  is  also  soluble,  in  the  LOW-field,  or  strong  coupling, 
sliding  regime. 

V.  The  Threshold;  a  Surviving  Challenge 

When  the  threshold  field  for  conduction  is  approached,  the 
incommensurate  chain  becomes  strongly  distorted  at  all  length 
scales.  As  a  result,  the  threshold  combines  the  all  the  richness 
of  critical  phenomena  (such  as  singular  responses  with  universal 
exponents)  with  the  additional  challenge  of  being  a  dynamic 

process,  far  from  equilibrium. 

The  ideas  of  the  renormalization  group  have  been  used 


20 


successfully  to  construct  a  transformation  which  has  a  fixed 
point  describing  the  depinning  transition  in  incommensurate 
dynamics9.  This  transformation  establishes  the  full  universality 
class  of  this  transition.  The  work  is  particularly  novel  and 
will,  it  is  hoped,  provide  a  basis  for  a  theory  of  the  threshold 
of  CDW  dynamics.  This  extension,  to  the  asymmetric,  dynamic  case 
of  the  conduction  threshold,  is  the  subject  of  ongoing  work. 


L.  Sneddon,  A.J.  Kassman,  S.  Liu  manuscript  in 
preparation . 


